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Non-Markovianity induced by a single-photon packet in a one-dimensional waveguide 
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The concept of non-Markovianity (NM) in quantum dynamics is still an open debate. Under¬ 
standing how to generate and measure NM in specific models may aid in this quest. In quantum 
optics, an engineered electromagnetic environment coupled to a single atom can induce NM. The 
most common scenario of structured electromagnetic environment is an optical cavity, composed by 
a pair of mirrors. Here, we show how to generate and measure NM on a two-level system coupled to 
a one-dimensional waveguide with no mirrors required. The origin of the non-Markovian behavior 
lies in the initial state of the field, prepared as a single-photon packet. We analyze how NM depends 
on two experimentally controllable parameters, namely, the linewidth of the packet and its central 
frequency. We relate the presence of NM to a 7r-phase shift between incoming and emitted fields. 
We also show how the two output channels of the waveguide provide distinct signatures of NM, both 
experimentally accessible. 

PACS numbers: 03.65.Yz, 03.67.-a 


I. INTRODUCTION 

The Markov approximation has been widely used to 
describe the dynamics of a quantum system coupled to 
its environment [T]. In this approach, the quantum dy¬ 
namics can be represented by a completely positive and 
trace-preserving (CPTP) dynamical map with a genera¬ 
tor in the Lindblad form - a quantum markovian master 
equation [5]. However, the Markov approximation can be 
violated when the coupling between the system and its 
environment is strong or when there is not a clear sepa¬ 
ration between the typical timescales associated with the 
system and the environment. Such a scenario has been 
reached experimentally, for example, in the contexts of 
quantum biology [3] and condensed matter physics [4], 
where the non-Mar kovian effects become relevant. 

Although the concept of non-Markovianity is well- 
established in the classical case, the definition of non- 
Markovianity in the context of quantum open systems is 
still subject of debate in the scientific community. For ex¬ 
ample, the definition proposed in Ref. [5] is based on the 
notion of non-Markovianity for a classical stochastic pro¬ 
cesses, namely, a quantum evolution is non-Markovian 
if it cannot be described by a divisible completely posi¬ 
tive trace preserving (CPTP) map. On the other hand, 
the non-Markovianity has also been characterized by a 
temporary back-flow of information from the environ¬ 
ment to the system 0 0. Other studies suggest that 
non-Markovianity can also be measured by the rate of 
change of the volume of accessible states of an open sys¬ 
tem [5], the amount of classical information extracted 
by the environment 0 , the appearance of negative deco¬ 
herence rates [TOj, and the smallest amount of isotropic 
noise that must be added to a quantum dynamics in or¬ 
der to describe it through a memoryless master equation 
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m- Besides its importance from a fundamental perspec¬ 
tive, the non-Markovianity can also be useful in applica¬ 
tions involving quantum communication [12] . quantum 
metrology [13] , quantum correlation generation just 
to name a few. 

The need to understand and to exploit non- 
Markovianity creates an interest in generating and ma¬ 
nipulating non-Markovian quantum evolutions. For ex¬ 
ample, Bi-Heng Liu et al. HD reported an all-optical 
experiment where non-Markovian aspects of the dynam¬ 
ics of a qubit - represented by the polarization degree 
of freedom of a photon - can be controlled by manipu¬ 
lating the initial state of its environment - represented 
by the frequency degrees of freedom. Another way to 
control the non-Markovianity of the dynamics of a qubit 
coupled to its environment is to manipulate the interac¬ 
tion between the qubit and an additional auxiliary qubit 
m, as demonstrated experimentally in a nuclear mag¬ 
netic resonance setup m- In the context of cavity quan¬ 
tum electrodynamics, the dynamics of a two-level atom 
may display non-Markovian effects due to the confine¬ 
ment of the electromagnetic field inside a cavity. This 
happens because the photon emitted by the atom can 
be reflected by the mirrors of the cavity and then reab¬ 
sorbed by the atom. Recently, T. Tufarelli and collabo¬ 
rators [18] showed that the dynamics of an atom coupled 
to a one-dimensional half-cavity - a semi-infinite waveg¬ 
uide with a perfect mirror at one end - can also provide 
a non-Markovian behavior. In this case, the presence 
of the mirror is also responsible for the origin of non- 
Markovianity. 

Here, we show that non-Markovianity in the two-level 
system (TLS) dynamics is induced by a single-photon 
packet in a one-dimensional waveguide. From a funda¬ 
mental perspective, there is novelty in the physical origin 
of non-Markovianity in our system. It is generated by a 
single-photon pulse, requiring no mirrors at the bound¬ 
aries of the TLS. We also provide the means to access 
the non-Markovianity of the TLS by making measures 
only on the field. From the applications perspective, our 
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approach allows for control of non-Markovianity by ad¬ 
justing the input state of the field, such as its central 
frequency and linewidth, requiring no change in the sys¬ 
tem parameters. 

This paper is organized as follows. In Sec. we in¬ 
troduce the model and solve for the equations of motion. 
Both the TLS and the field are studied. In Sec. m we 
explore the measures of non-Markovianity that are suit¬ 
able for our model. In Sec. |IV[ we present the signatures 
of non-Markovianity as they appear both in the TLS (A) 
and in the field (B) dynamics, both qualitatively and 
quantitatively. 

II. MODEL 

The model we consider (see Fig. consists of a 
one-dimensional (ID) infinite waveguide coupled to a 
two-level system (TLS) [T9H22] . whose ground and ex¬ 
cited states are denoted respectively by \g) and |e). 
The Hamiltonian of the TLS is Hs = hvscr+a-, where 
(T+ = \e){g\ = (7^1 and i/s is the transition frequency. 
The TLS is coupled to a ID electromagnetic environment 
composed by a continuum of frequency modes iz. Those 
modes are separated into two channels, a and b. Chan¬ 
nel a (b) describes the forward (backward) propagating 
modes Ui, {b^) with momentum = +\k^\ {kn = — \ki,\). 
The momentum is determined by the dispersion relation 
V = c\kv \ = ckjy. The Hamiltonian of the field is iLfieid = 
^ system-environment in¬ 

teraction Hamiltonian in the rotating-wave approxima¬ 
tion (RWA) is 

OO 

[^+ - h.c] , 

( 1 ) 

where rg is the TLS position, g^ is the coupling strength 
between the TLS and the mode with frequency v, and 
h.c. is the Hermitian conjugate. The RWA is valid in the 
weak coupling regime, that is, gv ^ ^s- In the following, 
we assume the reference frame where rg = 0. The total 
Hamiltonian is 

H = iLs+infield Tilint. (2) 

Since the number of excitations is conserved, the quan¬ 
tum state of the total system is written in the subspace 
of zero and one excitation, 

OO 

IV’(t)) = co(t)| 5 , 0 )-F?/’s(t)|e,0)-f ^<()^“)(t)4|5,0) 

OO 

+ ( 3 ) 

where |0) is the vacuum state, ^/’s(t) is the excited state 
probability amplitude, is the probability am¬ 

plitude of a photon with frequency v in the channel a 



Figure 1: (Color online) One-dimensional waveguide coupled 
to a TLS placed at rs = 0. The ground and excited states 
of the TLS are denoted respectively by \g) and |e). us is 
the transition frequency and Fid is the spontaneous emis¬ 
sion rate of the TLS. (6„) denotes the forward (backward) 
propagating modes of the electromagnetic environment with 
momentum k^, = +\k^\ {ki, = —\k^\). The TLS, initially in 
the ground state, scatters a single-photon packet described by 
the initial state 0). is the typical time duration of 

the pulse. The forward (backward) propagating packet after 
the scattering is measured by detector a{b) located at position 
rd{-rd)- 

(&), and Co(t) is the probability amplitude correspond¬ 
ing to no excitation. The dynamics of the total sys¬ 
tem is obtained by solving the Schrbdinger equation, 
ihdt\'4^{t)) = iL|'0(t)), that leads to the differential equa¬ 
tions 

V 

^(a).(fc)(^) ^ (5) 

co(i) = 0. (6) 

The last equation implies that co(t) is time independent. 

We investigate the case in which the TLS is initially 
in the ground state, V's(O) = 0- We also assume that the 
initial state of the field, in the real-space representation 
(/)(“)’(^)(r, t) = 'Yliv [22], is given by 

,),(“)(r,0) = iV 0(-r) 

(r,0) = 0, (7) 

where vg is the central frequency of the photon, A is the 
packet linewidth and iV is a normalization constant. The 
typical time duration of the pulse is (see Fig. [^. 
This initial state can be generated, for example, by the 
spontaneous emission of an auxiliary TLS, initially in the 
excited state 


A. Two-level System Dynamics 

The reduced dynamics of the TLS can be obtained 
in the Wigner-Weisskopf approximation [T]. In this ap¬ 
proach, we assume the continuum limit f dupiD, 
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where piD is the density of modes, and ^ gs- By sub¬ 
stituting the solution of Eq. ^ into Eq. Q , the dynamics 
of the excited state probability amplitude tps{t) becomes 


V’s(t) = - ( ^ + ws ) ipsit) 


( 8 ) 




where Fid = ‘inggpiB is the spontaneous emission rate of 
the TLS in the ID waveguide. Using the initial conditions 
from Eqs.Q, with N = y/Urpu^, and ^s(O) = 0, the 
solution of Eq. ® can be written as 


ipsit) = - 


FidA -('fiE. 


+ivs ) t 


--iSl 


- 1 




( 9 ) 

where Sl = — us is the detuning. Note that, in the 

monochromatic limit A ^ Fid, the excited state popu¬ 
lation tends to |'!/'s(t)P < 2 A/Fid. The behavior of the 
field in this limit is discussed in the next subsection. 

The reduced density operator of the system, ps{t) = 
Trfieid {\ip{t)){tp{t)\), in the {|e), |g)} basis, is 


Ps{t) = 




( 10 ) 


The reduced dynamics is given by the Master Equation 

m 


J^Psit) = -^S{t)[a+a_,ps{t)] 


( 11 ) 


+ r(t) (^a-ps{t)a+ - ^ {cr+cr_,ps(t)}^ , 


where 


the solution of Eq. ^ and applying the Wigner-Weisskopf 
approximation we find that 


—ct,0) 

+ \/nriuPiDQ{r)Q (^t - tps (t - 


(14) 


and 


(r, t) = \/7rFiDPiD0(-'r)0 (t -f V'S (^ + ^) ■ 

The first term in Eq( 14), namely (r — ct, 0), describes 
the free propagation of the initial wave packet. The sec¬ 
ond one describes the amplitude of a photon emitted by 
the TLS into channel a. The Heaviside theta functions 
guarantee that the emitted photon in channel a is de¬ 
tected at r > 0, after the photon arrived, t > rjc, at 
the detector at position r. Note that, for vanishing TLS- 
waveguide coupling Fid —>■ 0, the outgoing field is exactly 
the input field propagating in the direction defined by 
channel a at speed c, without changing the packet shape. 
For the counter-propagating packet, only the 

amplitude of the emitted photon is present, since the 
initial photon state contains no amplitude in channel b, 
i.e., + ct, 0) = 0. The Heaviside theta functions in 

Eq.(|T^ guarantee that a photon emitted into channel b 
is detected at r < 0, after it has arrived, t > |r|/c, at 
the detector at position — |r|. Eqs.(14) and (15) recover 
a well-known effect of light-matter interaction in one¬ 
dimensional environments, namely, total light reflection 
pn PSHU] . This is due to a destructive interference be¬ 
tween emitted and incoming fields [JSIESHSI]. It happens 
at resonance, 6^ = 0, and in the quasi-monochromatic 
regime, A Fid, for which [35] 


> 0, t) 


{5i:,=0, Ac^Fid} 


0 . 


(16) 


(H) 

represents the decay rate of the system, being Re(.) the 
real part, and 

(ID) 

represents the Lamb-shift, being Im(.) the imaginary 
part. Both quantities depend on time. 


B. Field Dynamics in the Real Space 
Repr esentat ion 


In the ideal monochromatic regime, no light will propa¬ 
gate through channel a after having interacted with the 
TLS, being totally reflected into channel b. Note that 
the origin of such destructive interference is the tt— phase 
shift between the incoming and the emitted packets. The 
field emitted into channel b also carries this 7r-phase 
shift with respect to the incoming packet, given that 
0 (fc)(|. < 0 ,t > |r|/c) cx ips(.t — kl/c) oc —(()i“((r — ct, 0 ). 

III. MEASURE OF NON-MARKOVIANITY 

In recent years different measures have been proposed 
in order to characterize and quantify the non-Markovian 
aspects of the dynamics of an open quantum system [ 51 - 
unida. When the dynamics of the system is described 
by a time-local master equation 


The dynamics of the field is studied in the real-space 
representation, as defined in the previous section. From 


( 17 ) 
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where 

C{t) = (18) 

with 




Tr[L,(t)]=0 e Tr 4(t)L,(t) 


= <5r 


(19) 


is a time-dependent generator of the dynamics in canon¬ 
ical form [3 [To], all measures agree that the quantum 
dynamics is Markovian if the time-dependent decay rates 
are always non-negative, that is, 7 i(t) > 0. In this case, 
the generator (181 is in Lindblad form [5] for each fixed 
t > 0. On the other hand, if the decay rates assume nega¬ 
tive values, the dynamics of the system is indivisible, that 
is, the quantum dynamics cannot be described by a se¬ 
quence of infinitesimal CPTP quantum evolutions 0171. 
According to the measure proposed by Rivas, Huelga and 
Plenio (RHP) 0, a quantum dynamics is non-Markovian 
if it corresponds to an indivisible CPTP quantum dynam¬ 
ics. 

Since the dynamics of the system of interest, as given 
by Eq. (11), is described by a time-local master equa¬ 


tion, the temporary appearance of a negative decay rate 
will be used as a witness for non-Markovianity. In other 
words, the TLS dynamics is non-Markovian if, and only 
if, r(t) < 0 for some t > 0. According to Ref. 0, the 
amount of non-Markovianity is given by 


■^= f{t)dt, 
Jo 


( 20 ) 


with f{t) = maxjO, — r(t)} > 0, where r(t) is the decay 
rate given by Eq. m. For the quantum dynamics de¬ 
scribed by Eq. ([TT]), the negativity of the decay rate can 
also be used as a signature of the non-Markovian behavior 
according to the measures proposed by Lorenzo et al. 0 , 
Hall et al. [ini, and Wolf et al. [U. Regarding the mea¬ 
sure proposed by Breuer et al. 0 , the non-Markovianity 
is not necessarily guaranteed when r(t) < 0. 

In our case, the non-Markovianity can be characterized 
by the dynamics of the excited state population ll/'s(t)p. 
To show this, note that Eq.(I2) can be rewritten as 


m = - 


dli’siOl" 




dt 


( 21 ) 


Thus, the decay rate becomes negative, r{t) < 0, when 
there is an increase in the excited state population, 
d\il^s{t)\'^/dt > 0 . 


IV. RESULTS AND DISCUSSIONS 


In this section we show the signatures of non- 
Markovianity evidenced both in the TLS dynamics (A) 


and in the field dynamics (B). In subsection A, we first 
present the existence of oscillations in the population of 
the excited state. We also show that these oscillations im¬ 
ply negative values for T(t), revealing non-Markovianity. 
The physical origin and meaning of such oscillations is 
explained. Finally, we quantify the amount of non- 
Markovianity in terms of the detuning 6^ and the packet 
linewidth A. In subsection B, we explore the means to 
directly access both r(t) and the TLS excited state pop¬ 
ulation experimentally, by making measurements only on 
the field. 


A. Signatures of Non-Markovianity in the 
Two-Level System Dynamics 

The dynamics of the excited state population of the 
TLS, 1^5(t)p, is plotted in Fi g0 (a). The analytical ex¬ 
pression is obtained from Eq.(|^. We set Tid = 1 as our 
frequency unit in all plots. The packet linewidth is kept 
at A = 0.1. The curves correspond io 5 l =0 (black 
dotted line), Sl = 2 (red dashed line) and Sl = 10 (blue 
solid line). The frequency of the population oscillation 
increases with the detuning, so as the number of peaks. 

We emphasize that the TLS is not in steady state for 
times TiDt ~ 8. In Figj^(b) we plot the same quantity as 
in figure (a), but for a time interval that is extended until 
Fid! = 50, keeping A = 0.1. The two curves correspond 
to Sl =0 (black dotted) and = 10 (blue solid line). 
Note that the latter has been multiplied by a factor 100 
so it fits the same vertical axis scale. In both cases, it is 
clear that the TLS is left in the ground state after the 
photon packet amplitude at the TLS position vanishes. 
This is illustrated in the inset. The time of interaction 
between the photon packet and the TLS scales as A“^. 

Fig0(c) shows T{t) for the same time interval as in fig¬ 
ure (a)) i.e., until Fid^ = 8. Again, the same detunings 
have been chosen, namely, Sl = 0 (black dotted line), 
Sl = 2 (red dashed line) and Sl = tt) (blue solid line). As 
discussed in section HI, at any given time t for which the 
excited state population is increasing, d\tps{t)\‘^/dt > 0, 
the decay rate becomes negative, r(t) < 0. It is particu¬ 
larly noticeable at higher detunings, e.g., at Sl = 10. In 
Sec. HI, we stated that F(t) < 0 implies the presence of 
non-Markovianity for the master equation of our model. 
Hence, population oscillation, or more precisely popula¬ 
tion increase, is a signature of non-Markovianity in the 
TLS dynamics. 

The physical origin of the oscillations in the TLS pop¬ 
ulations shall be explained. Firstly, it must be made 
clear that these are not Rabi oscillations |32|, neither 
classical [3T| nor quantum [33]. It is not classical be¬ 
cause the input is a quantum pulse, containing a single 
photon. It is not quantum Rabi oscillations either, since 
the photon emitted by the TLS should be otherwise re¬ 
absorbed after having been reflected by the mirrors of 
a cavity. In the present scenario, the emitted photon is 
never reabsorbed. The excited state amplitude dynam- 
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Figure 2: (Color online) (a) Excited state population of the 
TLS, \'tps{t)\'^, as a function of time for A = 0.1 and = 0 
(black dotted line), Sl = 2 (red dashed line) and 5l = 10 
(blue solid line). The population exhibits an oscillatory be¬ 
havior as the detuning increases, (b) Dynamics of the excited 
state population \'tps{t)\^ for a time interval that is extended 
until Fid! = 50, keeping A = 0.1. The black dotted curve 
corresponds to =0 and the blue solid line corresponds to 
Sl = 10. The curve associated with Sl = 10 has been multi¬ 
plied by a factor 100 so it fits the same vertical axis scale. The 
inset shows that for ct = 50c the photon packet has passed 
completely through the TLS. Thus, after Tiot = 50, the TLS 
reaches the steady state, where the excited state population 
becomes zero, (c) Time-dependent decay rate r{t) as a func¬ 
tion of time for A = 0.1 and the same values of Sl used in 
(a). The negative values of the decay rate, r(t) < 0, indi¬ 
cate that the TLS dynamics is non-Markovian. Note that 
r{t) <0 when there is an increase of the excited state pop¬ 
ulation, d\'ips{t)\'^/dt > 0. We set Fid = 1 as our frequency 
unit in all plots. 

ics, dtipsit), exclusively depends on the component of 
the photon packet that precedes the TLS, ct,0), 

not on > 0,t) neither < 0,t), which prop¬ 

agate outwards the waveguide. The remarkable feature 
of those oscillations is their dependence on the detun¬ 
ing. The number of peaks is proportional to 
This fact indicates that this effect is phase dependent, 
suggesting that it can be written in terms of an in¬ 
terference between two amplitudes, Ai and A 2 . Inter¬ 


ference between those two amplitudes is quantified by 
I = \Ai + A 2 P — (l^ip + |^ 2 p)- In Sec|n| we analyzed 
the interference effect in the field, between incoming and 
emitted photon amplitudes. On the contrary, we now 
seek amplitudes describing different states of the TLS it¬ 
self. For that purpose, Eq.(|^ can be recasted as dt^Js = 
^ 1 +^ 2 , where we define Ai = —(riD/2) V^s(t) and A 2 = 
—-\/(rA)/2 ^(“^(t) exp(—idit). Note that (j^^\ct,0) = 0 
due to the chosen initial condition. A rotating frame has 
been applyied, 'i/'s(^) = 6xp(—wst) and = 

exp (—At/2). We rewrite the increase of the excited state 
population, > 0, as 9t|Aip > 0. It follows that 

= -^^(2|Ai|2 + A;Ai+A2AJ) 

= -^(2|Ai|2+I). (22) 

The interference term I in Eq.( |22[ ) can be interpreted 
as a quantum interference between the amplitudes of 
the TLS to have already been excited (Ai) and to have 
not been excited yet (A 2 ) hence being in the ground 
state. In the latter case, it is the incoming field that 
carries the excitation. If no interference takes place, 
1 = 0, Eq.([^ reduces to spontaneous emission be¬ 
havior, 9t|AiE = —Tid |Aip. No increase of popula¬ 
tion is expected in this case. Only destructive interfer¬ 
ence, I < 0, can induce increase of the excited state 
population. We now analyze the dynamics of the in¬ 
terference term, X{t). Close to resonance for a quasi- 
monochromatic packet, <5^ ^ A ^ Tid, the dynamics of 
interference is X{t) « —2 |Ai(t)||A 2 (t)| Re[exp(— 
Therefore, in that regime, interference is periodically de¬ 
structive with frequency 5l, as qualitatively expected. 
Destructive interference becomes even more evident pre¬ 
cisely at resonance, Sl = 0 (A ^ Fid), for which 
ips{t) oc —(/^“^(—ct,0). In that case, Ai = |Ai| = —A 2 
and X{t) « —2 |Ai(t)||A 2 (t)| < 0. Such tt— phase 
shift between these two amplitudes is the same as the 
one responsible for zero transmission of the light field, 
(j)A){r > 0, t) « 0, as discussed in SecM Even in the res¬ 
onant finite-width packet regime A ;$riD there is still a 
peak in the population curve, as shown in Fig[^(b). In 
that regime, we have that Ai = |Ai| and A 2 = —IA 2 I, 
so destructive interference is maintained. However, after 
a given time t, the excited state population amplitude 
|Ai(t)| becomes bigger than the incoming-packet ampli¬ 
tude |A 2 (t)| so the system excitation declines, dt\Ax\^ = 
—Fid I All (|Ai| — IA 2 I) < 0. For the large-detuning 
quasi-monochromatic regime, Jj;, ^ Fid A, we have 
thatl(t) « —|Ai(t)||A 2 (t)| Re[exp(—i(i5Lt+ 7r))]. Again, 
interference periodically becomes destructive with fre¬ 
quency Sl^ so the peaks in Figj^(a) are spaced by a time 
interval of nearly l-nJS l- This last result is valid for both 
the low and the high detuning regimes. We have, there¬ 
fore, evidenced the physical origin of the successive in¬ 
creases in population, connecting them to tt— phase shifts 
that appear periodically. 

We now quantify the amount of non-Markovianity M 
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in terms of the detuning 5l and the packet linewidth 
A. Figj^(a) gives the detuning dependency of the non- 
Markovianity, N{5l). The curves are obtained by the 
use of Eq.(pO|). Three linewidth values have been chosen, 
A = 0.1 (black dotted line), A = 2 (red dashed line) and 
A = 10 (blue solid line). 


time is left for the population of the system to oscillate 
during the field passage. However, an even higher de¬ 
tuning will eventually rise non-Markovianity by inducing 
faster oscillations than the time spent by the packet to 
go through the TLS. 



Figure 3: (Color online) (a) Non-Markovianity A/” as a func¬ 
tion of detuning Sl for A = 0.1 (black dotted line), A = 2 
(red dashed line), and A = 10 (blue solid line). Note that, 
for Sl ^ A the non-Markovianity rises almost linearly. This 
is because the excited state population oscillates rapidly un¬ 
der this condition, leading to an increase of the time integral 
Eq.( |20[ |. (b) Non-Markovianity A/” as a function of linewidth 
A for Sl = 0 (black dotted line), Sl = 2 (red dashed line), 
and Sl = 10 (blue solid line). For Sl 0, non-Markovianity 
is maximum at A = 1. As we adopt Fid = 1, the maximum 
occurs when A = Fid. In this regime, light-TLS interaction 
is enhanced, inducing stronger oscillations in the excited state 
population. 


Three regions can be identified in the plot, corre¬ 
sponding to distinct behaviors. At very small detunings 
(5 l «C A, non-Markovianity Af is nearly constant. This 
is due to the tiny deviation from the ideal case of perfect 
resonance. The system dynamics is almost unchanged 
in that case. At intermediate detunings ~ A, a de¬ 
crease of Af becomes apparent. The amount of time dur¬ 
ing which there is increase in population gets reduced 
as detuning grows. That is the reason for the reduc¬ 
tion of non-Markovianity. At large detunings Sl ^ A, 
non-Markovianity rises almost linearly. As it has been 
discussed, large detunings induce fast population oscilla¬ 
tions. The successive peaks in the excited state popula¬ 
tion yield large Af when integrated in time. This third 
region where Af grows is not illustrated in the plot for the 
largest linewidth, A/Tid = 10. In such case, the photon 
packet is so short in time (and space), ~ A“^, that no 


Figj^(b) shows the linewidth dependency of the non- 
Markovianity, Af{A), for the detuning values of <5 l = 0 
(black dotted line), Sl = 2 (red dashed line) and <5^ = 10 
(blue solid line). The most prominent feature is the peak 
at A = 1, at non-zero detunings. The size of the peak 
increases with detuning. It is worth reminding that the 
unit here is the TLS natural linewidth Fid. So, A = Fid 
precisely corresponds to the mode matching condition 
between the linewidth of the incoming photon and that 
of the TLS. Therefore, light-TLS interaction is enhanced, 
inducing stronger oscillations in the excited state popu¬ 
lation. The peaks are asymmetric. In the large linewidth 
limit, non-Markovianity vanishes, whereas in the small 
linewidth limit, non-Markovianity is finite. Note that 
Af does not vanish in the A/Fid —t 0 limit even at 
resonance, because of the unavoidable initial increase in 
Ifsit 0^)P) starting at zero. Interestingly, the peak 
is absent for Sl — 0. In that case, oscillations in the ex¬ 
cited state population are minimized, as previously dis¬ 
cussed. The monotonic decrease of Af with increasing A 
at zero detuning is reasonable, since non-Markovianity 
is inversely proportional to the excited state population, 
Eq.(|2Tl). 


Our results are strongly based on the initial condition 
we have chosen. In Ref. [15], a different initial con¬ 
dition is addressed on a very similar context, namely, 
the TLS in the excited state and the field in the vac¬ 
uum state inside a ID waveguide. In that case, non- 
Markovianity takes place due to a single mirror put close 
to the atom. Non-Markovianity is then shown to vary 
with respect to the atom-mirror distance. In our model, 
on the other hand, the control on the degree of non- 
Markovianity comes from the choice of the shape and 
the frequency of the photon packet. This allows for dy¬ 
namical control without the need to change the system 
parameters, yielding a possible experimental advantage. 
We finish this section by calling attention to the impor¬ 
tance of a genuine single-photon packet for inducing non- 
Markovianity. Usual experiments with ID waveguides 
performed at the single-photon level employ very low- 
excitation laser pulses [Mj • Weak laser pulses can always 
be modeled by time-dependent Hamiltonians [32j . which 
in turn induce unitary dynamics. However, unitary dy¬ 
namics are Markovian, since negative decay rates do not 
appear, as shown in Eqs.(I8) and (20 1 . Therefore, for 
our system of interest, an authentic single-photon packet 
in the input field is necessary for generating finite non- 
Markovianity. 
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B. Signatures of Non-Markovianity in the Field 
Dynamics 


So far we have investigated the TLS dynamics. From 
an experimental perspective, accessing the system dy¬ 
namics can only be done by measuring the emitted field. 
In ID electromagnetic environments, however, there can 
exist interference between the incoming and the emitted 
fields, as previously mentioned. So we need to carefully 
specify what kind of information is encoded in the field 
propagating towards each of the two detectors placed out¬ 
side the waveguide. The intensity of a photodetection 
signal is /a.b(t) = where e is the vac¬ 

uum field constant and is the position of the detector. 
This relation is in agreement with Glauber’s photodetec¬ 
tion theory [55] . 

Firstly, we analyze the signal from the backward- 
propagating packet (channel b), detected at position 
Tfi < 0. Eq.(15) shows that provides direct ac¬ 


cess to the excited state population of the TLS, Ib{t) oc 
kd|/c)|^. Therefore, the TLS be¬ 
havior would be directly measured as shown in Figj^(a). 
The blinking of the TLS in the monitor of detector b 
would indicate the existence of non-Markovianity, man¬ 
ifested by successive increases and decreases of popula¬ 
tion. 

The signal from the forward-propagating packet (chan¬ 
nel a), however, exhibits richer phenomena. Eq.(14) evi¬ 


dences the interference between the incoming packet and 
the one emitted by the TLS. Sections |ll| and |IV| present 
detailed analysis on the consequences of that interfer¬ 
ence. Being more specific, it has been shown that the 
TT—phase shift is responsible for both the total reflec¬ 
tion of light and the increase of the excited state popula¬ 
tion, that evidences non-Markovianity. Hence, destruc¬ 
tive interference induces asymmetry between backward 
and forward propagating fields on the one hand and non- 
Markovianity on the other hand. It is natural, then, to 
look for a signature of non-Markovianity on the ratio be¬ 
tween and Indeed, it can be 

shown that 


ID 


Re 


_(/)W(-|rd|,t)_ 


= F ( t - ) , (23) 


where the term F(t), on the right side, is defined in 
Eq.([l^. Non-Markovianity is exactly quantified by its 
negativity. It is worth mentioning that S(t — |rd|/c) = 
-2^S where S{t) is the 

time-dependent Lamb-shift also defined in Eq.([I^. Re¬ 
markably, Eq.(|2^ provides the means to directly obtain 
the curve shown in Figj^(c) in an experiment, by exclu¬ 
sively measuring the outgoing electromagnetic fields in 
the appropriate manner. 


V. CONCLUSIONS 


We have shown that a single-photon wave packet gen¬ 
erates non-Markovianity in the dynamics of a two-level 
system coupled to a one-dimensional waveguide. Both 
the TLS and the field dynamics have been obtained. We 
derived an effective master equation for the TLS open 
dynamics. The form of the master equation allowed us 
to measure non-Markovianity in terms of the negativ¬ 
ity of the time-dependent decay rate. The decay rate is 
shown to be negative if and only if the time derivative 
of the excited state population is positive. Oscillations 
in the excited state population are explained in terms 
of periodic 7r-phase shifts that appear between incom¬ 
ing and emitted fields. Non-Markovianity Af is quanti¬ 
fied in terms of packet linewidth and detuning. Af in¬ 
creases almost linearly at large detunings and presents a 
peak at the mode matching condition A = Fid. Both 
linewidth and detuning are experimentally controllable 
parameters, yielding a practical advantage. The signa¬ 
tures of non-Markovianity can be accessed by measuring 
the light at both ends of the waveguide. The intensity 
of the backward propagating field is shown to be propor¬ 
tional to the excited state population of the TLS. The 
asymmetry between forward and backward propagating 
fields reveals non-Markovianity, providing information on 
the time-dependent decay rate. 

We highlight two perspectives offered by this work. 
From an experimental point of view, measuring spectra 
is generally easier than measuring time resolved quan¬ 
tities. Therefore, relating non-Markovianity of the TLS 
dynamics to the spectra of the light field in both chan¬ 
nels of the waveguide can be of high interest. From a 
theoretical point of view, it is relevant to understand the 
influence of other quantum states of the field on the non- 
Markovianity of the TLS. In particular, if the TLS is 
initially excited and a photon packet is prepared in the 
waveguide, we could relate stimulated emission to non- 
Markovianity. In that case, analyzing Breuer’s measure 
of non-Markovianity would be suitable. Another interest¬ 
ing problem is to analyze the TLS dynamics after a two- 
photon packet is prepared inside the waveguide. Tempo¬ 
ral correlations between the two photons could generate 
unexpected NM behaviors. 
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